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ABSTRACT 

Recent observations have found the ratio of Ha to free-free radio continuum to 
be surprisingly high in the diffuse ionized ISM (the so-cahed WIM), corresponding to 
an electron temperature of only ~3000 K. Such low temperatures were unexpected in 
gas that was presumed to be photoionized. We consider a 3-component model for the 
observed diffuse emission, consisting of a mix of (1) photoionized gas, (2) gas that is 
recombining and cooling, and (3) cool H I gas. This model can successfully reproduce 
the observed intensities of free-free continuum. Ha, and collisionally-excited lines such as 
[N Il]6583. To reproduce the low observed value of free-free to Ha, the PAH abundance 
in the photoionized regions must be lowered by a facto r ~3, and ^20% o f the d iffuse 



Ha must be reflected from dust grains, as suggested by lWood &: Reynolds! (119991 ) . 



Subject headings: atomic processes, ISM 



1. INTRODUCTION 

Low density ionized regions, often referred to as the warm ionized medium (WIM), account for 
Vo or more of the ionized interstellar hydrogen in the Galaxy. The WIM occupies a substan- 
tial volume fraction and is a major component of the interstellar medium (ISM) at ^1 kpc aboy e 



the galactic disk plane in both the M i lky Wav (IMcKee 
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Although it is now generally 
believed t hat only niassive 0-type stars have enough ionizing power to account for the ionization in 
the WIM ([Revnoldsl 12004 ) . the detailed ionization mechanism is still not fully understood. Struc- 
tures such as superbubbles and chimneys have been invoked to explain how H-ionizing photons 
could travel from the stellar sources to ionize the WIM. 

Observations of optical diagnostic lines such as Ha, H/3, He IA5876, [N II] A6548, 6583, 5755, 
[S II]A6716, 6731, [O I]A6300, [O II]A3727, and [O III]A5007, as well as infrared lines ([Ne II]12.81 fim 
and [Ne III] 15.55 /im) h ave been made in the la s t two decades, aiming to understand the physical 



conditions in the WIM ( Greenawalt et al. 



1997 



Randlll997l . l200d : iHaffner et al.lll999l : lOtte et al 



2002; iMiller &: Veilleuxl l2003l : iMadsen et al.l |2006| ) . Although the observed line ratios vary from 
region to region, optical line ratios such as JN II]A6548/Ha appear to indicate that the WIM 
material has temperature T w 8000 it 2000 K. iHill et al.l (120081 ) estimate the characteristic electron 



density in the WIM to be Ue ~ 0.07cni~^, with a volume filhng fraction of ~30%. Optical line 
ratios such as [O III]A5007/Ha indicate that compared with H II regions powered by early 0-type 
stars, the WIM seems to have fewer ions p resent tha.t require ioiiization energies greater than 23 eV, 
with n(He"'")/nHe ranging from 0.3 to 0.6 (IReynoldslllOSSl : lRandll997l : iHaffner et al.lll999l : iRevnolds 
2OO4I : iMadsen et a l. 2006: iHaffner et al.ll2009l ). 



Recently, iDobler et al.1 (J2009|, hereafter DDF09), using the WMAP five-year data, deduced 
a surprisingly low temperature (~3000 K) for the WIM from the measured ratio of the thermal 
bremsstrahlung (free-free) emission to the associated Ha emission. Similar res ults had previously 
been reported based on the free-free to Ha ratio in the WMAP one-year data (JDavies et al.1 120061 ) 
and three- year data (JDobler fc: Finkbeineri l2008l ) . The low value came as a surprise, because the 
diffuse ionized gas had been expected to have T w 8000K, based on both observed optical line 
ratios and the predictions of steady-state photoionization models. 

Previous discussions of the temperature and ionization of the WIM have generally assumed 
steady-state conditions. Here we instead argue that much of the WIM is not in equilibrium, and 
that it is possible to understand the observed line ratios and weak free-free emission if part of the 
WIM is assumed to be gas in the process of cooling and recombining after removal of a photoionizing 
source, as a consequence of either stellar evolution or changes in the opacity along the sightline 
between the parcel of gas and the O star providing Lyman continuum photons. We propose that 
observations of free-free and line emission from the diffuse ISM at high latitudes typically sum 
over multiple components: (1) gas currently being photoionized, (2) gas that was photoionized 
in the past but which is currently cooling and recombining, and (3) neutral gas (the cold neutral 
medium and warm neutral medium) . Dust grains mixed with the gas also contribute scattered light 
that includes emission lines from HII regions in the disk. For plausible choices of parameters, we 
find that this model can account for the low ratio of free- free to Ha found by DDF09, while also 
reproducing observed optical line ratios, such as [NII]6583/Ha. 

The structure of this paper is as follows. In Section [2] we describe the model used to simulate 
cooling and recombination in the WIM. The simulation results are presented in Section [3l where 
we explore various factors which affect this process. In section H] we propose a model using three 
ISM components to explain the observed line ratios and free-free emission. We discuss the main 
result and our models in Section [5l followed by a brief summary in Section [6l 



2. Model Description 

We model the evolution of temperature, ionization, and emission from cooling, recombining gas. 
As initial conditions we take the ionization and temperature to be consistent with photoionization 
by a distant OB association, but at t = we assume that the radiation with hu > 13.6 eV is suddenly 
turned off. The remaining ionization and heating sources are cosmic rays and photoelectric emission 
by the diffuse interstellar radiation field with hv < 13.6 eV. We follow the ionization of 11 elements: 



H, He, C, N, O, Ne, Mg, Si, S, Ar, and Fe, which includes all the elements that normally have gas 
phase abundance relative to H above 10~^, and includes all important coolants. 



2.1. Physics in the Model 



The evolution is assumed to take place at constant volume, with H nucleon density tt-h = const. 
The abundance XA,r = n{A~^^)/n}i of ion A^"^ evolves according to 



dx y^ 7- 
dt 

dxAfl 
dt 



CA,r-lXA,r-^l + 



(rr) , (dr) \ , (qr) 



A , / (»*»■) I {dr)\ , (or) 



. (rr) , (de)\ , (or) 

He I a^l 1 + ttyl 1 I + "-H"^ 1 



XA,r 

XA,i - Ca,oxa,o 



XA,r+l 
for r > 1 



(1) 
(2) 



where (^A,r is the probability per unit time of ionization A'^^ — )• A'^^^^ + e due to either photoion- 
ization, cosmic rays, or secondary electrons; a^^ is the rate coefficient for radiative recombination; 

ttjJJ is the rate coefficient for dielectronic recombination; and aJ^'J is the effective rate coefficient 
for recombination on dust grains. The present model includes only ions with r < 2, and gas tem- 
peratures T < 10 ^ K. We assu me case B recombination for H and He. The recombination rate for 
H is taken to be (|Drainell201ll ) 



aB = 2.54xl0-13j.-0.8163-0.02081nT4^^3g-l ^ 

where T4 = T/IO'^K, and for recombination of He II— )-He I we take 

as (He) = 2.72 x iQ-i^r^-o.Tsg ^^^3 g-i _ 



(3) 



(4) 



Fo r othe r elements, radiative r ecombination rat e s are evaluated by subroutine rrfit.f from 
Verneii (119991 ). whic h uses rates from IPeauienot et al.l (119911) for ions of C, N, O , and Ne (refitted 
with the formula of lVerner fc Ferlandl ( 19961 )). and from IShull fc van Steenberg ( 1982 ) for ions o f 



Mg, Si, S, Ar, a nd Fe. Rate coefficients a^'^^> fo r Mg and Si were from lNussbaumer &: Storey 



for S were from IShull &: van Steenbergi (| 19821 ). and for Fe were from 



Arnaud Sz Ravmond 




Rate coefficients a^'l for recombination on grains are taken from IWeingartner &: Draind ( 2001al ) 



The electron density Ue = nn "^a Z^r '''XA,r, and the free particle density n = 11^ + nn ^^ ^^ xa,: 



For isochoric evolution, the temperature T evolves according to 



dT 
H 



A 



(3/2)n/c 



^1 dn 
n dt 



(5) 



where (F — A) is the net rate of change of thermal kinetic energy per unit volume due to heating 
and cooling processes (including ionization and recombination), and dn/dt is the net rate of change 
of the free particle density due to ionization and recombination processes. The heating rate per 



volume r includes heat deposition by cosmic ray ionization and by photoelectrons emitted from 
atoms, ions, and grains. The cooling rate per volume A includes kinetic energy removed from the 
gas by inelastic collisions and by rec ombination of ions and electrons. The mean kinetic energy per 
recombining electron is taken to be (|Draindl201ll . eq. 27.23) 



{E„) ^ [0.684 - 0.0416 In r4] kT 



(6) 



Radiative cooling processes consist of free-free emission, and line emission following collisional 
excitation of various atoms and ions. Collisional deexcitation is included, although it is unimportant 
at the densities considered here. 

Grain-related processes are included in our simulation. From their study of the "spinning 
dust" emission from the WIM, DDF09 concluded that the PAH abundance in the WIM is a factor 
~3 l ower than in the H I. Because the PAHs account for a large fraction of the photoelectric heat- 
ing ( Bakes &: Tielenslll994l : IWeingartner &: Draind 2001br). and are also thought to dominate the 
grain- assisted recombination ( Weingartner Sz Draind l2001al ). depletion of the PAHs will affect the 
heating and recombination. To explore this, we multiply the rates for dust photoelectric heating 
and grain- assisted recombination by a factor g, where g = 1 gives the rates estimated for normal 
grain abundances in the diffu se ISM for photoelectric heating and grain- assisted recombination 
(jWeingartner &: Draind l2001al ). Our standard model for the WIM assumes g = 1/3, which pre- 
sumably results primarily from reduced abundances of the small grains that account for most of 
the grain surface area. We will explore the sensitivity to the reduction factor g by also performing 
simulations with g = 1 (no reduction in PAH abundance) and g = 0.1, a factor of 10 suppression of 
grain photoele ctric heating and grain-a,ssisted recombination. Grain photoelectric heating is calcu- 
lated following I Weingartner Sz Draind ( 200 Ibl) . and the g = 1 rates for grain- assisted recombination 
are taken from 



Weingartner &: Draind (J2001a 



We calculate collisional excitation and resulting radiative cooling by ions of C, N, O, Ne, Si, 
S, Ar, and Fe. We include a total of 159 lines, but the 26 lines of listed in Table [T] account for more 
than 95% of the radiative cooling at each point in the thermal evolution of our models. Sources of 
the collisional rate coefficients for species used in our calculations are listed in Table [2j 



Table 1: Principal Cooling Lines (wavelengths in vacuo) 



[Cll]157.7/im 
[C II]2328A 
[C II]2326A 
[N Il]205.3^m 
[Nll]121.8^m 
[N II] 6585 A 
[N II]6550A 



[Ol]145.5/im 
[O I]63.19/xm 
[O I] 6302 A 
[O II]3730A 
[O II]3727A 
[Nell]12.81^m 
[Si Il]34.81^m 



[S II] 6733 A 
[S II]6718A 
[S II] 4070 A 
[S IIl]33.48/im 
[Slll]18.71/im 
[S III] 9533 A 
[S III] 9071 A 



[Fe Il]25.99/im 
[Fe Il]5.340/^m 
[FeII]1.644//m 
[Fell]1.321//m 
[Fell]1.257/um 



Table 2: Sources for Collisional Rate Coefficients 



Ion 
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H" 
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2pO 
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'Pj- 


-'B2 


01 


'Pj- 


-'Pj' 


01 
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-ID2 


on 


400 

■^3/2 


-2D° 


Nell 


2po 

^3/2 


2pO 


Sill 


2po 

-^1/2 


2pO 

- ^3/2 


SII 


4qo 
*3/2 


-^D^j 


SIII 


^Pj- 


-'Pj' 


SIII 


'Pj- 


-'B2 


Fell 


'Bj- 


-'Bj, 


Fell 


'Bj- 


-'Fj' 


Fell 


'Bj- 


-4Dp 



Taval (2008) 
Hudson fc Bell f2005 ) 
Hudson fc Bell ('2005 1 
Peauignot (1996) 
Pequignot (1996') 
Taval (2007) 
Griffin et al. (2001 1 
Bautista et al. (2009) 
Taval &: Zatsarinnv (2010) 
Taval fc Gupta (1999) 
Taval fc GuDta (1999 ) 
Ramsbottom et al. (2007 ) 
Ramsbottom et al. (2007) 
Ramsbottom et al. (2007) 



Barinovs et al. (2005) 



Abrahamsson et al. (2007) 



Barinovs et al. (2005) 



For the diffuse interstellar radiation field we use the estimate of iMathis et al.l (|l983l ). Pho- 



toionization rates for this radiation field were taken from 



using photoionization cross sections from IVerner fc YakovlevI (119951 ) and IVerner et al.l (jl996l ) . 



Draind (2011, Table 13.1), calculated 



We include secondary ionization and heating by cosmic rays following iDalgarno &: McCrav 



(119721). Charge exc hange between H and O is included in our model, with rate coefficients from 
Standi et alT(|l999l ). 



2.2. Free-Pree vs. Ha 



We calculate the free-free emission at 41 GHz, Ij^(41GHz), to compare with observations 
(DDF09): 

> = 5.44 X 10-^^ gsT^^-^e-'"'/''^nl erg cm^ s"^ Hz"! sr'^ (7) 



where Ue is the electron num ber density, and the Gaunt factor (|HummeiJll988l ) is accurately ap- 
proximated by (JDraindl201ll . eq. 10.9) 



gs ~ In 



Vs. 



exp 5.96 - ^— IniugT^:'-'') + e 



where vg = z^/GHz. The Hq emission rate is: 

Ma = 2.82 X 10-26r-°-^^2-°-°^^'"(^*)nen(H+)erg cm^ s'^ sr'^ 



(8) 



(9) 
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Fig. 1. The ratio of jrV(41GHz)/(JHa/hi^) as a function of gas temperature T for H ionization fraction xjj = 1 (solid 

curve), 10~^ (dash curves), and lO"** (dot curves). For each ajjj, upper curves employ gas phase elemental abundances for 
-F* = to simulate WIM, lower curves employ gas phase elemental abundances for Fi, = 0.8 to simulate CNM. Observed nearly 
all-sky averaged free-free/Ha ratio ~ 0.85kjy sr~^ R~^ from WMAP (DDF09) is indicated by the horizontal dash-dot line, 
corresponding to gas temperature ~3000 K for fully ionized gas. 



corresponding to an effective rate coefficient 



47rJHc 



_ - T T„ ^ ,r,-l3r7.~0.942--0.0311n(T4) 

"Ho = — ; = 1-17 X 10 T. cm s 

hu 



3„-l 



(10) 



The ratio of free-free to Ha depends on both the temperature and the fraction of the electrons 
contributed by H"*": 



V-m 



> (41 GHz) 

JRa/hu 



0.0465 



Up 



-10.442+0.031 In r4 



n(H+)' 



In 



exp j 3.913 + -^ IutA + e 



kJy sr ^ 
R 
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Free- free emission comes from all ions, therefore in low-ionization gas we have free- free emission 
from ions such as C"*" + e~ , whereas Ha comes only from H"*" + e~. Therefore, if H becomes 
almost neutral, tp is sensitive to the gas-phase abundance of elements that can be ionized by 
hu < 13.6 eV photons. Figure [1] shows ip{T) as a function of gas temperature from 50 — 10^ K, 
for three different H ionization fractions xh = n(H"'")/nH: xh = 1 (solid curve), 10~^ (dash curves) 
and 10~^ (dot curves), and for two gas phase elemental abundances: F.^ = 0, representing WIM, 
and F^ = 0.8, representing CNM. The horizontal dash-dot line indicates the observed free-free/Ha 
ratio ~ 0.085kJy sr^^ R~^ determined from the WMAP 5 year data (DDF09). The observed ratio 
corresponds to gas temperature ~3000 K for fully ionized hydrogen. 

In addition to emission from the WIM, there will be some Ha and free-free emission from 
H I clouds. The H I gas in the interstellar medium is found at a wide range of temperatures. The 



majority is in the CNM phase, at T ss 10^ K, with ne/n(H+) ss 2 (if CcR/n{R) ~ 2 x lO^^'^cm^s"^ 
then ~50% of the free electrons are from C^, S^, and other metal ions.) From this gas we expect 
ip ~ 0.017[ne/n(H"'")] kJy sr~^ R~^ ~ 0.034kJy sr~^ R~^. In addition, a substantial fraction of the 
H I is in the WNM phase, with T w 10^ K; this gas will have ip ^ 0.047 kJysr^^ R^^. Overall, we 
estimate that the Ha and free-free emission from the H I phase will have ^hi ~ 0.04. 



2.3. Initial Conditions 



Gas-phase elemental abundances are based on the recent s t udy b y I Jenkina (|2009l ). For our 
standard model, we use abundances from the model of I JenkinsI (|2009l ) with F^, = 0, representing 
a relatively low level of depletion (see Section [5] for a discussion of the choice of F^). For N e 
and Ar (not covered in Jenkins' study), we assume solar abundances from lAsplund et al.l (J2009l ). 
He/H=0.1 is also assumed. In addition, we construct a model with F^, = 0.25 (reduced gas phase 
abundance for elements that deplete) to exp lore the sensitivity to coolant abundances. Carbon 
requires special attention, since recent work (ISofia &: Parvathil |2010| ) indicates that the oscillator 
strength of C II]2325A is larger than previously estimated, implying that gas phase C abundances 
estimated from measurements of C II] 2325 A could be lowered by a factor ~2. For this reason, in 
addition to the standard model and the F^, = 0.25 model, we consider the "Reduced C" model in 
which C has an abundance two thirds of its F-*, = value and other elements all have their F^, = 
values. 

For our standard model we take nn = 0.5 cm~^ and initial temperature Ti = 8000 K. Under 
these conditions the thermal pressure of the photoionized state, p/k ~ 2.15nHr ~ 9 x 10'^ cm~^ K, 
is comparable to, although somewhat higher than, current estimates for pressures in the diffuse 
ISM. In our stan dard model, the initial ionization fractions (IIF) of all the elements (Table [3]) 
are adopted from ISembach et al.l (|200Cl ). where we use the va lue in their standard model with 
parameter Xedge = 0.1 for fully ionized gas (ISembach et al.ll200Cl . Table 5), corresponding to a fairly 
soft ionization field. Two other IIFs (Table [3]) considered include values from the Orion nebula 



([Baldwin et al.lll99ll ). which corresponds to a hard radiation field , and a set of observed values 



based on various previou s studies of the WIM (JHaffner et al.lll999l : iRevnoldsl l2004l : iMadsen et al 



20061 : iHaffner et al.ll2009l ). The effects of initial gas density and temperature are studied as well. 
Table |4] summarizes all the models in our simulations. 

The cosmic ray primary ion ization rate for an H atom, CcR) is controversial. Direct mea- 
surement of the cosmic ray flux (jWang et al.l |2002| ) at £■ < 1 GeV is affected by the solar wind. 
Webber &: YushakI (119831) att empted to correct for solar wind modulation. Smoothly extrapolating 



the 



Webber &: YushakI (|l983l ) spectrum to lower energy su ggests Ccr, ^ 1-3 x 10 ^'^ s ^. However 



studies of the ionization conditions insid e molecular clouds ( Black &: van Dishoecklll99ll : lLepplll992l : 
McCall et al.ll2003l : llndriolo et al.ll2007l ) indicate primary ionization rates CcR ~ 0.3 — 3 x 10~^^ s~^, 



and we will consider values of CcR within this range. We adopt CcR = 1 x 10 s for the standard 



model, and will then explore varying (^CR from 2 x 10 



-17 



-1 to 5 X 10" 
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0.01 
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Fig. 2. — EfTect of different density and cosmic ray ionization rate. All cases have g — 1/3 and F* = 0. Top 
panel: The temperature evolution for nn = 0.5 cm^'^ (solid curves) and nu — 2.5 cm~'^ (dashed curves), and different 
cosmic ray ionization rate (en- Bottom panel: The ratio of cumulative free- free emission at 41 GHz to cumulative Ha 
emission for models with different ("cr for nn ~ 0.5 cm~^. Values of xu = n(H"^)/nH are given at several points along 
the curves by different tick marks. Note that when Ccr/'^h > 6 x 10~^® cm'^s"^, the final temperature Tf > IQ^K, 
and ^{nut = 10^ cm'^^ yr) > O.OQkJy sr"^ R'^ 

3. Results 



The temperature as a function of time for our standard model is shown in the top panel of 
Figure [21 Initially the gas cools rapidly, reaches a minimum temperature ~100K within 0.4 Myr, 
and then warms up to 400 K, approaching its asymptotic state in several Myr. The rise in temper- 
ature at nut ^ 2 X 10^cm~^yr occurs because the [CII]158/im coolin g declines as free electrons 
recombine. Similar late time reheating has been seen previously (e.g., lDraindll978l ). Define the 
cumulative emission ratio as: 



^{t) 



L{jna/hv)dt' 



(12) 



^ starts from ~0.17 kJy sr^^ R^^, corresponding to the initial temperature Tj = 8000 K, then 
drops as the gas cools, to ~0.08 for 10^ cm^'^ yr < nnt ^ lO^cm^'^ yr as shown in the bottom panel 
of Figure [2j This low value of ^ is consistent with the observed value of ~ 0.085kJy sr^^ R^^ 



(DDF09). 

Altering the parameters in the standard model will change the cooling history, and we explore 
the influence of various factors. Based on their impact, we divide the influential factors into two 
classes: major factors, which includes CcRj ^Hj abundance of metal elements and grain depletion, 
and minor factors: the initial temperature and IIF. 



3.1. Major factors 

Figure [2] shows the sensitivity to CcR and tt-h- For fixed tt-h, the final asymptotic temperature 
(Tf) increases with increasing cosmic ray ionization rate, while for fixed CcRj Tf decreases with 
increasing nn- At a higher density the final temperature and free- free to Ha emission ratio become 
less sensitive to CcR- For ^CR ranging from 2 x 10^"'^'^ s^^ to 5 x 10^^^ s^"*^, Tf increases from 
~240 K to ~3800 K when nu = 0.5 cm"^^ i-,^^ Qj^jy fj.^^ ;^20 K to 190 K when nu = 2.5 cm'^ (as 
shown in top panel of Figure [2]). As CcR is varied from 2 x 10^^^ to 5 x 10~^^ s^^, the emission 
ratio ^(nnt = 10^cm~^yr) changes about 20% for nn = 0.5 cm^^ (as shown in bottom panel of 
Figure [2]), but only about 1% for nu = 2.5 cm~^. For each density, there is a critical value of CcR, 
above which the gas remains warm, with asymptotic temperature T > 3000K, and ^ > 0.10. 

Elemental abundances affect the cooling process. Figure [3] shows the effect of varying the 
element abundances. In general, reducing gas phase abundances will reduce the cooling, increase 
the final temperature Tf, and consequently will raise ^. As we discussed in Section [2. H evidence 
from analyzing spinning dust emission in the WIM (DDF09) suggests that the small polycyclic 
aromatic hydrocarbons (PAHs) may be underabundant in the WIM relative to the general ISM. 
We explore the effect of different grain reduction factors, with g = 1/3 as our standard model. As 
shown in Figure [3l the final value of ^ is quite sensitive to the value of the PAH reduction factor g. 
When g is reduced, the photoelectric heating rate drops; since the grain-assisted recombination rate 
also drops, more electrons and ions will be present in the gas phase, leading to increased cooling 
through collisionally-excited lines. The combined result of the two factors is that the ionized gas 
recombines more slowly, cools faster, and reaches lower values of ^ when g is small. 



3.2. Minor factors 

The initial temperature and ionization fractions also affect the cooling process. Quantities like 
T and xh for models with different initial T and ionization will evolve differently at the initial 
stage of the simulations, but will all converge to the same asymptotic steady-state values. On the 
other hand, the cumulative ratio of free- free to Hq, ^, will differ because it involves integration 
over time. We study the effect of different Tj by running a model with T, = 10"^ K (instead of 
8000 K), and we study the e ffect of different IIF by running models with IIF2 and IIF3. IIF3 



assumes Orion Nebula values ([Baldwin et al.lll99ll ). formed by a much harder ionization field than 



10 




n^t(10''cm"V) 

Fig. 3. — EfTect of varying gas-phase abundances (F* = 0.25 vs. F* = 0), PAH abundance {g — 1 vs. g — 1/3), 
and C abundance (Reduced C vs. F* =0). Top panel shows the temperature as a function of time for different 
models (indicated by model name), and bottom panel show the ratio of cumulative free- free emission at 41 GHz to 
cumulative Ha emission as a function of time. 

the one that f orms the WIM, as we discussed above in Sect i on [H IIF2 is a set of values based on 



the hterature (jMadsen et al.l |2006J; iHaffner et al. 



199S 



20091 : lRevnoldsll2004l ). fahing between IIF3 



and our standard model. These three models cover a large range of IIF. As shown in Figure HJ the 
results are insensitive to changes in the different initial ionization conditions. 



4. A Three Component Model for the Diffuse Emission 



In this section we try to compose a model to simultaneously explain the observed low free-free 
to Ha emission ratio (corresponding to T ~ 3000 K), and line ratios such as [N II]A6583/Ha and 
[N II]A5755/[N II]A6583 (indicating T > 8000 K) in the diffuse emission from the WIM. We note that 
the observed low free- free to Ha emission ratio (DDF09) comes from measurements which integrate 
over a large fraction of the high-latitude sky, while the measurement of various other diagnostic 
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Fig. 4. — EfTect of different initial conditions (initial temperature and ionization fractions) on the evolution of 
temperature T and the cumulative ratio \& of 41 GHz free-free to Ha. 

lines from the WIM, like [N II]A6583 or [N ll]A5755, are limited to relatively high intensity regions. 
In our model, the observed diffuse emission comes from three components: 



1. "hofphotoionized gas with T ~ lO^K, nearly fully ionized, producing free-free emission, 
recombination radiation, and collisionally-excited lines such as [N II]A6583 and [N II]A5755; 

2. cooling and recombining gas, which starts from the "hot" photoionized phase and cools and 
recombines to give Ha and free-free emission as well as a small amount of collisionally-excited 
lines, such as [N II]A6583; 

3. neutral H I gas, partially in the "cold neutral medium" (CNM), with T « lO^K, and partially 
in the "warm neutral medium" (WNM), with T « 10^ K; together these emit a small amount 
of Ha and free-free emission due to cosmic ray ionization, but negligible metal line emission 
in the optical. 
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In addition to the actual emission from the H I gas, there will also be reflected light - both Ha 
and metal lines such as [N II]A6583 - that was originally emitted elsewhere, mainly regular H II 
regions, and then reflected by dust grains present in th e H I g as. Let f^^ be the fraction of the 



observed Ha that is scattered light. IWood Sz Revnoldsl (|l999l ) estimated that f^^ « 0.05 — 0.20 



P(refl) 



at high galactic latitudes. In our calculation, we assume the reflection coefficients of [N II]A6583, 
[N II]A5755 and [S II]A6716 are the same as Ha, because the wavelengths are close to Ha. 

The observed intensities are weighted averages over the three components. Here we ask what 
weighting factors are needed to reproduce the observed low ratio of free-free/Ha, as well as other 
line ratios. 

The relative contributions of the three components can be determined from the observed emis- 
sion ratios. Among all the diagnostic line ratios in the WIM, [N II]A6583/Ha and [S II]A6716/Ha 
are the two best-stu died cases, while other lines have been studied only in a few select directions 
(JHaffner et al.ll2009l ). These two line ratios depend on temperature and ionization fraction of N and 
S in the gas. The high second ionization potential of N (29.60 e V) protects it frorn being doubly 
ionized in the WIM, which has a generally soft ionization field ( Madsen et al.ll2006l ). This factor 
along with the similar fir st ionization potent i als of N and H (14.5 3 and 13.60 eV) make N"'"/N close 
to unity in the H II gas (JMadsen et al.ll2006l : iHaffner et al.lll999l ). On the other hand, some of the 
S in H II regions will be doubly ionized due to its low second ionization potential (23.33 eV), which 
leads to uncertainty in the predicted [S II]A6716/Ha, reducing the utility of [S II]A6716/Ha as a 
constraint. Nevertheless, we have calculated the integrated value of [S II]A6716/Ha for the cooling 
gas (see Tabled. 

Let f^° and f^^ be, respectively, the fraction of the observed diffuse Ha emitted from the 

be the fraction of the observed Ha that is actually 

(refl) 



f(refl) 



photoionized gas and H I clouds, and let /jj^ 

reflected from dust in H I. Then f^°° ™^ = 1 — f^° — f^^ — /^q"'' is the fractional contribution of 
the cooling material. The emission from this three component model depends on four parameters: 
and the temperature Thot of the photoionized component. For an adopted 



„(hot) „(HI) „(refl) 
JUa ' JKa ' ■'Ho 



value of /jj^ , a physical solution must have < /jj^ < 1 



< fSL'^ < 1 



J. (refl) „ ^ „ (cooling) ^ . 



/, 



(refl) 
H« ' 



and 



fua ■ ^^ ^^^^ ^^^ observed free-free/Ha and [N II]A6583/Ha emission ratios as 



two constraints to solve for the relative fraction of the three components. 
The instantaneous [N II]A6583/Ha line ratio is: 



HT) 



J[NII]A6583 _ -^2 4 2^0.495+0.040 In T4 -2.204/T4 ^(N^)/^(H+; 
JHa ^ ■ ' 7.41x10-5 



(13) 



where we use the collision strength Q('^Pq,^ D2) 

(|2005l ). The cumulative ratio of [Nil] 6583 to Ha for cooling gas is 



Q3Q3j.0.0528+0.009inT4 frpm iHudson k Belli 



Ht)^ 



,fo i[NII]6583 dt' _ /o Ht') JHa dt' 



/o jUa dt' 



f JRadt' 



(14) 
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Fig. 5. — Ha fractions for the three gas components based on two constraints (1,^(41 GHz) /Hq — 0.85kJysr ^R ^ 
and [NH]A6583/Ha=0.4), as a function of hot (photoionized) gas temperature Thot- The equations are solved under 
these assumptions: f^^f = 0.20, f^J^ = 0.014, n{}i+)/n^ = 0.5 in the HI (corresponding to Ccr ~ 1 x 10"^^s~^ in 



CNM), and He+/He=0.3 in the hot ionized gas, as in HFl. The fact that f^^ ' + f^ 



f(HI) 



/jj^° is coincidental. 



Cooling gas in our standard model gives $ ~ 0.06 (see Table [3D. The WNM and CNM contribute 
negligible [NII]A6583. 



We take the av e rage value of 



Revnolds et al. 



2001 



Madsen et al. 



N II] 6583 /Ha = 0.4, based on observations (JRevnolds 



2004 



20061 ). Then we must have 



p(hot) ^(rcfl) 



0.4 = (/C^ + C O0(rhot) + (1 - /^ 



(hot) 
a 



.(HI) 
J Ha 



fSf)'^ 



(15) 



Note that in this equation, we assume the reflected component has the same [N II]A6583/Ha ratio 
as the hot ionized component, as discussed previously. 



The ratio of all-sky free-free to Hq emission is /j^(41 GHz)//(Ha) 
(DDF09). Thus we must have 



0.085 kJy sr-i R- 



0.085 kJysr^^R"^ 



/^^"^^(rhot) + (1 



„(hot) 
JUa 



.(HI) 
J Ha 



.(refl)x,^, , .(HI) , 
fna )^ + /Ha V'HI 



(16) 



For trial values of Thot and /^^ > 0, we can use the observed /([NII]6583)//(Ha) and /jy(41GHz)//(Ha) 
to determine f^° and /^^ , by solving the two linear equations ifTS]) and (fT6]) . The coefficients 
ip(T]^ot) and (/'(Thot) are obtained from eq. (fTT|) and (fTBj) . while $ ^ 0.08 and $ ^ 0.06 are ob- 
tained from the appropriate simulations as their asymptotic values (see Table [Has well as Figure [2]). 
Based on the discussion in §2.2| we take V'Hi ~ 0.04. Physical solutions must have < f^^ < 1, 



0</i^'^<l,and0<(l- 



.(hot) 
JUa 



.(HI) 
JHa 



■ /jj^ ). Figure [5] shows one example, where the reflected 



Ha fraction is set to 20%, ?7-(H"'")/ne = 0.5 in the H I phase and n(He^)/n(He) 
gas phase. The sensitivity to the flrst two parameters will be discussed below. 



0.3 in the hot 
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Fig. 6. — (a) Temperature Thot and (b) fractional Hq contributions from different ISM components, as func- 
tions of the reflected Ha fraction f^^ ', based on two constraints — 7^(41 GHz)/Ha = 0.85kJy sr~^R~^ and 
[NII]A6583/Ha = 0.4 — for fllf^ ranging from 5% - 30%. It is assumed that only f^J^ ^ 0.014 of the Ha is 
emitted by the CNM and WNM components (see text). Panel (a) also shows the dependence of the solution (Thot) 
on the value of n(H^)/ne in the HI gas. Our standard model assumes ^cr = 1 x 10~^^s~^, which corresponds to 
n(H^)/ne ~ 0.5 in the "cold" phase, as discussed in Section|4l 

In addition to the requirement that the fractions of all the three components have to be positive, 
which, for f^^ ^ 0.2, requires 8300< Thot < 9800 K (see Figure [5]), there is a second restriction 
coming from the fact that the total amount of Ha emission from the H I phase is determined by 
the cosmic ray ionization rate. Consider a uniform layer of cold neutral gas with a column density 
A^(HI), with cosmic ray ionization of H balanced by case B recombination as well as grain assisted 
recombination; the Ha intensity from this component at latitude b is 



hv 



1 3 X lO^o cm 



-2 



47r 



sin 6| 



-Ccr(1 + 4>cr) 



QHa/«B 
1 + agr/XeOCB 



(17) 



where Xe = n^/nH, (t>CK ^ 0.67 is the number of secondary ionization per primary ionization in 
neutral gas (IDalgarno fc McCravlll972l ). and A'^(HI) ^ 3x 10^*^ | sin6|^^ cm~^ ([Radhakrishnan et al 



1972 : iDickev et al.lll978l l. Thus 



hv 



0.01 

I sin6| 



CcR 



lO-i^s-i 



R 



(18) 



where the H I is taken to be a 1:1 mixture of CNM and WNM material 



a 



The distribution of Ha for \h\ > 10° has recently been measured bv iHill et al.l ( 20081 ). who 
found the full WIM to be fitted on average by lua/hv ^ (0.625±0.002)| sin b\~^ R. Comparing these 
two results, the fraction of Ha emitted by the H I should only be f^^ w 0.014((^cr/10~^^s~^). 



^We have taken (^cr ~ 0.67, QHc/as « 0.6, a^r/xeas ^ 10 for CNM with T ^ 100 K, n^ ~ 0.01, nn ^ 30 cm^ 



( Draindlioiil ): ^cr ~ 0.5, ana/as f» 0.5 and ag^/xeaB ^ 0.5 for T ~ 5000K WNM; and A^'c: 
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Fig. 7. — [NII]A5755/[N IIJA6583 averaged over components (solid curve) as a function of reflection fraction based 
on the solution in panel (b) of Figure [T] Dashed curve and dotted curves show [NII]A5755/[NII]A6583 associated 
with pure cooling and pure hot gas phase, giving low er and upper limits of the ratio on individual sightlines. The 
result agrees with observations in lMadsen et al.l (|2006r ). although comparison between the two is not straightforward 
(see discussion in Section |4|. 



In Figure [6] we take /jj^ = 0.014 and consider different values of j^^ . For each /jj^ , the 
Ha/free- free and [N II] /Ha constraints serve to determine /pj^ and Thotj as shown in Figure El 
For /jj^ ranging from 5% to 30%, panel (a) in Figure [6] gives T^ot ranging from 8000-15000 K. 
For /i'f^ = 0.2, we find Thot ~ 9100 K. 

One uncertainty when solving equations ()15p and (J16p is n(H"'')/ne in the H I, which affects 
■^Hi (equations [TT]) . Assuming ^CR ~ 10~^^s~^, Tcnm ~ 10^ K and ncNM ~ 30cm~^, about half of 
the electrons come from metal elements, so the free-free/Ha ratio for the H I will be ^hi ~ 0.04, as 
discussed in Section Yl% We explore the effect of this uncertainty on T^ot in the range 0.02 — 0.06, 
corresponding to (n(H+)/ne)coid = 1 and |, as shown in panel (a) of Figure [6l We find this 
uncertainty does not significantly affect the derived solutions. 

In addition to [N Il]A6583/IIa, the [N II]A5755/[N II]A6583 ratio is another indicator of the 
temperature in the ionized medium. While [N II]A6583 is one of the strongest lines in the WIM, 
comparable in intensity to Ha, [N II]A5755 is extremely hard to detect due to its low intensity. 



Never theless, several measurements of this line ratio have been reported (Table 6 in lMadsen et al. 



(J2006l )). The [N II]A5755/[N II]A6583 data have a large scatter, ranging from 0.023 (corresponding 
to a temperature over 12000 K) down to the observational upper limit 0.002 (corresponding to a 
temperature below 6000 K). In our WIM model the [N II]A5755/[N II]A6583 ratio is determined 
by the relative amounts of emission from the hot and cooling components, with the ratio for the 
cooling gas about 0.004 (almost model independent, as shown in Table Hj). The solid curve in 
Figure [7] shows the weighted [N II]A5755/[N II]A6583 ratio based on the solutions in the panel (b) 
of Figure E] for j^^ ranging from 5%-30%, under the assumption that the [N II]A5755/[N II]A6583 
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ratio of the reflected light is the same as its value in the hot ionized gas, as we discussed previously. 
Note that this result is a weighted average value over all the components in the WIM - measure- 
ments along individual sightlines will vary. The observed [N II]A5755/[N II]A6583 value along a 
specific sightline depends on the weight of the three components along this sightline, which may 
differ from the weights that give the observed all sky averaged free- free to Ha and [N II]A6583 to 
Ha ratio. The dashed curve and the dotted curve in Figure [7| show [N II]A5755/[N II]A6583 of 
the cooling gas and the hot gas, which could be interpreted as the lower and upper limit of the 
average [N II]A5755/[N II]A6583 value under extreme weighting cases. Nevertheless, the average 
weighted [N II]A5755/[N II]A6583 va lue ranges from 0.008 to 0.025, in general agreement with the 



observed values and upper limits in iMadsen et al.l (|2006l ) . Individual sightlines could have values 



of [N II]A5755/[N II]A6583 as low as 0.004, if the sightline is dominated by the cooling gas. 

In sum, if the cooling gas is described by our standard model, the temperature of the hot gas 
and the fractions of the three components could be solved for so that a weighted sum of fractions 
of the three components: 

1. reproduces the low observed free- free/Ha ratio; 

2. reproduces the observed [N II]A6583/Ha ratio 

3. has [N II]A5755/[N II]A6583 in agreement with observations; 

4. has a small fraction (~1.4%) of the Ha emitted from H Igas; 

5. includes a reflected component accounting for ~20% of the total Ha 

We emphasize again that the three-component model presented in this section is aimed to 
explain averaged observational results; individual sightlines will generally differ from the average 
because (1) the relative weights of its components differ from the average values, the temperature 
Thot and the ionization conditions in the H II component may vary, and (3) the cooling component 
along an individual sightline may not be an average over the entire cooling and recombination. 
There is evidence that the physical conditions in the WIM do v ary from sightline t o sightline, 



and even among different velocity components within one sightline (jMadsen et al.ll2006l ). Different 
regions in the WIM will have different densities, cosmic ray ionization rates, elemental abundances, 
and histories. The three component model presented here is highly idealized, but it appears to 
provide a physical framework that is consistent with the observations. 



5. Discussion 

When the hv > 13.6 eV starlight is cut off at the beginning of the simulation, the hot gas begins 
to cool and recombine. The temperature, particle density and gas pressure decrease. In general, 
the cooling time scale is short, as shown in Figure [2) for nn = 0.5 cm^^, T falls below 10'^ K within 
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0.3 Myr. If the photoionized gas was initially overpressured (relative to its surrounding), it would 
initially be expanding, resulting in adiabatic cooling after the photoionizing source turns off at 
t = 0. Conversely, if the photoionized gas was in pressure equilibrium with a confining medium at 
t = 0, it would begin to undergo compression as it cools and recombines at t > 0. However, if the 
nn = 0.5 cm~'^, photoionized region is ~5 pc or larger, the cooling time will be short compared to 
the sound crossing time, and the effect of adiabatic expansion or compression should be of secondary 
importance. In the present work we assume isochoric evolution, which appears to be a reasonable 
approximation during the cooling phase. 

In modeling the cooling gas, we assumed case B recombination. In the real WIM, some of the 
/iz^ > 13.6 eV recombination radiation might escape from the WIM and get absorbed by nearby 
cold or warm neutral gas (T < 1000 K); the resulting ionized H in these cooler media would 
recombine, emitting Ha and free-free emission, with a ratio of free-free emission to Ha appropriate 
to T < 1000 K, lower than the ratio for the cooling gas. This would have the effect of decreasing 
the ratio of free-free emission to Ha below the value calculated for the cooling gas model, allowing 
the observed low free-free/Ha ratio to be explained by a somewhat smaller value of /fj™°™^ i the 
fraction of the Ha contributed by cooling and recombining gas. However, we expect this "escape" 
of ionizing photons to be minimal: the mean free path of ~14eV photons in partially recombined 
H is short, and we expect the case B on-the-spot treatment to be a valid approximation. 

In our standard model, we choose -F^ = for the elemental abundances in the WIM. In 
Jenkins' model, the nonzero ini tial depletion (F^, = 0) is identified wit h the gas-phase abundances 



in a warm, low-density medium (jSpitzerJll985l : ISavage fc Sembachll 19961 ). and should therefore apply 



to the WIM. Also, as mentioned above, there is evidence (DDF09) showing that the PAHs in the 
WIM are underabundant by a factor of ~3. Our calculation shows that both factors are crucial 
for cooling the gas quickly enough and to a low enough temperature to be able to reproduce the 
observed low free-free/Ha ratio. As shown in Figure [3l if there is no grain depletion or there is 
appreciable depletion of coolants (such as for F^ = 0.25), the recombining gas stays relatively hot 
in the final steady state, and gives an unacceptable high integrated free-free/Ha ratio. 

Cosmic rays heat the gas directly, and they also raise the electron density n^. If rie is low, 
photoelectron emission causes dust grains and PAHs to become positively charged, reducing the 
photoelectric heating rate. Thus increased cosmic ray ionization, by lowering the charge state of the 
dust and PAHs, has the effect of increasing the dust photoelectric heating rate. For our cooling gas 
model to be able to reproduce the low free-free/Ha ratios that are observed, the ratio of cosmic ray 
primary ionization rate to gas density Ccr/^h should not exceed ~ 6x 10~^^ cm^ s~^. For Ccr/^-h ^ 
5 X 10~^^cm^s~^ (e.g, nu = 0.5 cm~^ and CcR = 5 x 10~^^s~^, with Ccr/^-h = 1 x 10~^^ cm^ s~^), 
the cosmic ray ionization maintains n(H"'")/nH > 0.1, and grain photoelectric heating can sustain 
the gas at T > 10^ K. Observation s of Hj" give estimat es of CcR ^ (0.5 — 3) x 10^^^ cm^ s~^, with 



an average of ~ 2 x 10 ^^ cm^ s ^ (jindriolo et al.l 120071 ). If CcR ~ 2 x 10 ^^ cm'^ s ^, then the gas 



density cannot be much smaller than tt-h = 0.5 cm ^ if the gas is to cool with ^ < 0.085, as required 
to explain the low observed value of free-free/H a. This high initial value of rig = 0.5 cm~^ seems 
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at first sight to be at odds with the study by iHih et al.1 (|2008l l. which concluded that the most 
probable value of n^ is only rig ~ 0.03 cm~^ in turbulent models of the WIM. Note, however, that 
in our standard model, the cooling gas ends up with xh < 0.03 and rig ~ 0.015cm~^. 



Because the observed line emission and free-free emission are weighted sums over 3 compo- 
nents, one of which (the cooling phase) itself has a range of temperature, there is no way to measure 
emission ratios associated with different phases unless the components could be separated. Tem- 
peratures derived from different line ratios need not agree. Moreover, even for a single line ratio, 
the contributions from different components may vary significantly from sightline to sightline, or 
with velocity on a single sightline. This may explain the large scatter in the physical conditions 
deduced from different line ratios. 

The present model envisages intermittent photoionization events in the diffuse high-latitude 
gas, followed by cooling and recombination. Suppose that the probability per unit time of a pho- 
toionization event is Tpi, with the ionizing radiation lasting a time thii before the photoionization 
switches off and the gas begins to cool and recombine. During the recombination phase, the number 
of Ha photons per recombining H"'" is ~0.5. The ratio of the Ha emission from the photoionized 
gas to the emission from the cooling gas is 



.(cooling) -1 ^ 0.5 

■I Ha P' 

Hence, 



~ ■'Ha c^tw 

rieana f (™°i'°g) 

•'Ha 

If f^f^ w 0.2, then /H^°*V/Hr"°^^ ~ 0.22/0.56 « 0.39, and the duration of the photoionization 
phase is only thii ~ 1 x 10^(0.5 cm^^/ne) yr. 

Given that O star lifetimes are ~ 3 x 10^ yr, such a short value of thii seems surprising. Longer 
values of thii can be obtained if /^'^ is larger: if f^^ = 0.3, then f\^° I fna' ~ 0.35/0.33 « 
1.1, raising thii to ~ 3 x 10^(0.5cm^'^/ne)yr, but this is still short compared to O star lifetimes. 
Longer values of thii can be obtained if rig is lowered, but if rig is much smaller than 0.5 cm~'^, 
grain photoelectric heating and cosmic ray heating prevent the gas from cooling to low enough 
temperature to be able to account for the observed low free-free/Ha ratio, if, as assumed, the 
cooling takes place at ~ constant density. 

O star lifetimes are not the only time scale on which photoionization can vary. The ionizing 
radiation for the WIM may be provided in large part by runaway O and B stars, with space 
velocities > lOOkms"^. Modulation of the photoionization rate at a given point may result from 
motion of the ionizing sources through a medium of varying opacity: a star moving at lOOkms"^ 
travels 10 pc in 10^ yr. 

Although O stars have been generally favored as the source of ionization for the WIM, it 
should be kept in mind that the actual source for WIM ionization has not been securely established. 
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According to the present analysis, the ionization phase, if it heats the plasma to ~ 10^ K, needs to 
be of relatively short duration in order to be able to explain the low observed ratio of free-free/Ha. 
Perhaps other transient processes (such as shock waves or soft X-rays due to supernova blast waves, 
infalling gas, etc.) contribute to heating and ionization of the WIM. 

A prediction of the current model is that the Ha emission should have a component with the 
same radial velocity profile as the 21 cm emission, resulting from cosmic ray ionization of the H I. 
However, because cosmic ray ionization is thought to account for < 2% of the Ha, it may not be 
easy to recognize this component. The reflected component of Ha, [N II], etc, should also correlate 
with A^(HI), but the apparent space velocities will differ from 21 cm radial velocities due to motion 
of the emitting H II relative to the reflecting dust grains. 



6. Summary 

The principal points of this paper are as follows: 

1. We simulated the cooling and recombination of initially photoionized gas in the WIM following 
removal of hi' > 13.6 eV radiation. The ratio of free-free emission to Ha, and emission line 
ratios such as [N II]A6583/Ha, were calculated, and various factors which influence the cooling 
were explored. The result strongly depends on the abundances of elements, gas density, 
cosmic ray ionization rate, and the abundance of very small grains and PAHs, but depends 
only weakly on the initial temperature and ionization fractions, as shown in Figure El [3] and 

a 

2. Based on these calculations, we propose a three component model — emission from hot 
ionized gas, cooling gas and neutral H I, plus reflected light ~ to explain multiple observations 
in the WIM. With plausible weighting factors for these components, the model simultaneously 
yields the low free-free to Ha ratio, which indicates a low temperature (~3000K, DDF09), 
and metal hue ratios, such as [N II]A6583/Ha and [N II]A5755/[N II]A6583, which indicate 



a high temperature (~ 10 K, iMadsen et al.l (120061 )). The reflected component is crucial 



there must be a fracti on (^20%) of the observed Ha coming from the reflected light, consistent 



with the estimate bv lWood &: Revnoldsl (|l999l ). to explain the observation 



3. For our model to successfully reproduce the low free-free/Ha ratio, some restrictions of the 
physical conditions in the WIM are required: 

(a) The ratio of cosmic ray primary ionization rate to gas density Ccr/^h should not exceed 
~ 5 X 10~^^cm^s~^ (Figure El), consistent with current observational estimates of CcR if 
rill ~ 0.5cm^^. 

(b) The gas phase element depletion parameter F^ < 0.25 (Figure ^ so that there are 
sufficient gas phase coolants to cool the gas in the pres e nce o f heating by cosmic rays 



and photoelectrons from grains, consistent with I Jenkins! ()2009l ). 
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(c) The abundance of ultrasmall grains (including PAHs) should be suppressed by a factor 
~3 relative to the abundances in the overall H I (Figure [3|), consistent with DDF09. 

4. Sightline-to-sightline variations of emission ratios are expected in our model. First, the pa- 
rameters which affect the cooling model, such as the gas density, elemental abundances, 
cosmic ray ionization rate and the depletion factor of small grains, may vary spatially within 
the WIM. Second, different sightlines will have different proportions of the three components 
in our model. Third, the cooling component is time-dependent, and individual sightlines may 
not present a full average over the entire cooling history. These effects will lead to the vari- 
ation of physical quantities deduced from observations on different sightlines, or for different 
velocity components on an individual sightline. 
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Table 3. Initial Ionization Fractions 



IIF = 1^ IIF = 2'' IIF : 



A 


UA/nii'^ 


I 


II 


III 


I 


II 


III 


I 


II 


III 


IV 


H 


1 


0.05 


0.95 




0.05 


0.95 







1 






He 


0.1 


0.68 


0.32 





0.5 


0.5 





0.04 


0.96 







C 


10-3.65 


0.01 


0.94 


0.05 


0.02 


0.83 


0.15 





0.13 


0.87 





N 


10-4.21 


0.05 


0.93 


0.02 


0.1 


0.8 


0.1 





0.13 


0.87 





O 


10-3.25 


0.05 


0.95 


0.00 


0.05 


0.95 


0.05 





0.43 


0.57 





Ne 


10-4.03 


0.11 


0.89 


0.00 


0.3 


0.69 


0.01 





0.91 


0.09 





Mg 


10-4.65 


0.01 


0.70 


0.29 


0.02 


0.63 


0.35 





0.05 


0.95 





Si 


10-4.61 


0.00 


0.87 


0.12 


0.01 


0.74 


0.25 





0.05 


0.83 


0.12 


S 


10-4.58 


0.00 


0.78 


0.22 


0.02 


0.68 


0.3 





0.03 


0.95 


0.02 


Ar 


10-5.56 


0.01 


0.76 


0.23 


0.15 


0.75 


0.1 





0.11 


0.89 





Fe 


10-5.41 


0.00 


0.34 


0.64 


0.01 


0.69 


0.3 





0.01 


0.24 


0.75 



''Standard Model in lSembach et al.1 l|200Cl 'l with Xedge = 0.1 

''Es timate based on lHaffner et al.l l|l999l '): lRevnoIdj l|2004l 'l : iMadsen et al.l l(200^ ): lHaffner et al.l 
20091 1: between IIFl and IIF3 

'^ Orion nebula values IIBaldwin et al.lll99ll ) 

"^Elemental abundance for standard Model (_F* = 0) 
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Table 4. Models for Recombining Gas 





Name'' 




(cm-3) 


Ccr" 
(10-i''s~i) 


F.d 


g° 


IIF* 




v]>h 


$i 


[NII]A5755 ■ 
[NII]A6583 


[SII]A6716 ^ 
Ha 




standard 




0.5 


1.0 





0.33 




8000 


0.081 


0.060 


0.0039 


0.15 




riH = 2.5 




2.5 


1.0 





0.33 




8000 


0.080 


0.063 


0.0039 


0.16 


"•H 


= 2.5 CCR = 


0.2 


2.5 


0.2 





0.33 




8000 


0.080 


0.063 


0.0039 


0.15 


riH 


= 2.5 CcR = 


2.0 


2.5 


2.0 





0.33 




8000 


0.080 


0.062 


0.0039 


0.15 


"H 


= 2.5 CcR = 


3.0 


2.5 


3.0 





0.33 




8000 


0.079 


0.062 


0.0039 


0.15 


nn 


= 2.5 CcR = 


5.0 


2.5 


5.0 





0.33 




8000 


0.079 


0.061 


0.0039 


0.15 




CcR = 0.2 




0.5 


0.2 





0.33 




8000 


0.081 


0.062 


0.0039 


0.15 




CcR = 2.0 




0.5 


2.0 





0.33 




8000 


0.083 


0.059 


0.0039 


0.15 




CcR = 3.0 




0.5 


3.0 





0.33 




8000 


0.088 


0.060 


0.0039 


0.15 




CcR = 5.0 




0.5 


5.0 





0.33 




8000 


0.111 


0.078 


0.0032 


0.42 




F* = 0.25 




0.5 


1.0 


0.25 


0.33 




8000 


0.099 


0.085 


0.0040 


0.14 




Reduced C 




0.5 


1.0 


0+C' 


0.33 




8000 


0.088 


0.065 


0.0039 


0.16 




9 = 1 




0.5 


1.0 





1.0 




8000 


0.109 


0.083 


0.0038 


0.27 




9 = 0.1 




0.5 


1.0 





0.1 




8000 


0.075 


0.056 


0.0039 


0.13 




IIF2 




0.5 


1.0 





0.33 


2 


8000 


0.083 


0.054 


0.0038 


0.14 




IIF3 




0.5 


1.0 





0.33 


3 


8000 


0.083 


0.028 


0.0028 


0.063 




Ti = 10* K 




0.5 


1.0 





0.33 


1 


10000 


0.086 


0.099 


0.0069 


0.23 



''Model name 

''H nucleon density 

'^Cosmic ray primary ionization rate. 

"^Depletion parameter from l|jenking|2009l ) 

''The grain reduction factor. The grain assisted recombination rate and photoelectric heating rate reduce to this fraction of their 
full values. 

'initial ionization fractions (see Tabled 

^Initial gas temperature. 

''Integrated ratio of free-free at 41 GHz to Ha, as in Equation 1111 for njjt = 10^cm~^yr. 

'Integrated ratio of [Nil] A6583 to Ha, as in Equation 1141 for nut = 10^cm~^yr. 

JIntegrated ratio of [Nil] A5755 to [Nil] A6583, for nut = lO'^cm-^yr. 

■^Integrated ratio of [SII] A6716 to Ha, for njit = lO'^cm-^yr. 

'Element abundance identical to the standard model (_F* = 0) except for C, for which we take the Fi, = value x(2/3), or 
nc/i^H = 1-4 X 10~^ 



